Ultrafast dynamics of neutral superexcited Oxygen: A direct measurement of the 
competition between autoionization and predissociation 
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Using ultrafast extreme ultraviolet pulses, we performed a direct measurement of the relaxation 
dynamics of neutral superexcited states corresponding to the nlag(c^l2~) Rydberg series of O2. An 
XUV attosecond pulse train was used to create a temporally localized Rydberg wavepacket and 
the ensuing electronic and nuclear dynamics were probed using a time-delayed femtosecond near- 
infrared pulse. We investigated the competing predissociation and autoionization mechanisms for 
superexcited molecules and found that autoionization is dominant for the low n Rydberg states. We 
measured an autoionization lifetime of r = 92it6 fs and r = 180 di 10 fs for (5s, 4d)ag and (6s, 5d)ag 
Rydberg state groups respectively. We determine that the disputed neutral dissociation lifetime for 
the V — {) vibrational level of the Rydberg series is Tq — 1100 =b 100 fs. 



A pervasive theme in ultrafast science is the charac- 
terization and control of energy distributions in elemen- 
tary molecular processes. Ultrashort light pulses are used 
to excite and probe electronic and nuclear wavepackets, 
whose evolution is fundamental to understanding many 
physical and chemical phenomena [1 . However, until 
recently, time-resolved studies of wavepacket dynamics 
using light pulses in the infrared (IR), visible, and ultra- 
violet (UV) regime had predominantly been limited to 
low- lying excited states and femtosecond timescales [2 . 

Advances in photon technologies, specifically in the 
field of laser high-harmonic generation (HHGjfS'.l , have 
opened new avenues in the time-resolved studies of molec- 
ular dynamics. The efficient energy regime for HHG lies 
within the extreme ultraviolet (XUV) range from 10 eV 
to 100 eV, allowing for the excitation of inner-shell elec- 
trons to form highly excited atomic and molecular states. 
Importantly, the attosecond to few-femtosecond nature 
of the HHG radiation forms a temporally localized ex- 
cited state wavepacket whose dynamics can be followed in 
real-time using a time-delayed probe [5-9 . Recently ex- 
periments have employed these attosecond XUV sources 
to study ultrafast fragmentation [10 and autoionization 
[ITJ [12] in highly excited molecular ions. 

Apart from excited molecular ions, neutral molecules 
existing far above the ionization potential form another 
important class of molecular systems where the appli- 
cation of ultrafast XUV sources could provide new in- 
sights into the fast relaxation mechanisms. These neu- 
tral superexcited molecular states are often found as 
intermediates in chemical reactions initiated by high- 
energy photons (e.g., solar radiation) and play an im- 
portant role in the chemistry of planetary atmospheres 
[T3] . The superexcited molecules quickly relax from their 
non-equilibrium state through multiple competing decay 
mechanisms including autoionization, fiuorescence, and 
dissociation into neutral fragments p!4HT7) . Synchrotron 
light sources have traditionally been employed for inves- 
tigating these processes. While these sources offer ex- 
cellent spectral resolution, they lack the time resolution 
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FIG. 1. An ultrashort XUV pulse excites a 2sau electron to 
form the nlag{c^ll~) Rydberg states of O2 (blue curves) con- 
verging to the c^S" state of (green curve). These curves 
represent quasi-bound superexcited molecular states that can 
only support two vibrational levels, = 0, 1. Since the states 
are imbedded in an electronic and dissociation continuum, au- 
toionization and neutral dissociation are the two prominent 
decay mechanisms for these highly excited molecules. 

required for the real-time study of ultrafast relaxation 
channels in superexcited molecules. In contrast, an at- 
tosecond pump-probe experiment can resolve dynamics 
on the few femtosecond timescale. 

In this work, we focus on the dynamics of the superex- 
cited oxygen molecule, specifically the nlcjg family of neu- 
tral Rydberg states converging to the c^S~ state of O2". 
These superexcited states are formed through the direct 
excitation of an inner-shell 2scFu electron to states that 
lie greater than 10 eV above the ionization potential of 
O2. The Rydberg state potential energy curves resem- 
ble the c^S~ curve and also exhibit quasibound vibra- 
tional levels, V = 0, 1. The potential energy curves for 
the (n = 4, h)s(7g{c^Ti~) states are shown in Fig. [l] They 
are of special interest because of their importance in the 
dynamics of ozone formation and decay in Earth's upper 
atmosphere [13]. This system presents a unique oppor- 
tunity to study the competition between neutral dissoci- 
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ation and autoionization and explore the non-adiabatic 
effects which dominate the dynamics of these superex- 
cited states. 

For the high n Rydberg levels, the dynamical evolu- 
tion is expected to mimic the c^S~ ion-core dissociation 
[18], which has been studied extensively through syn- 
chrotron experiments p^H2Q] and theoretical calculations 
[2TH23] . Tanaka et al [2T1 proposed tunneling dissoci- 
ation through the low potential barrier as a prominent 
decay mechanism and estimated dissociation lifetimes of 
To = 18 ns for z/ = and ri = 50 ps for u = 1. In the first 
experimental study of the dissociation dynamics, Evans 
et al. [19] measured dissociation lifetimes of tq = 270 ±30 
fs and Ti = 69±7 fs using pulsed field ionization at a syn- 
chrotron facility. Recently, Demekhin et al. [23] calcu- 
lated the diabatic and adiabatic potential energy curves 
of {c^T>~)02 and estimated the dissociation lifetimes of 
To = 12 ps and ri = 68 fs, which agrees with the neu- 
tral particle data taken by Hikosaka et al. [18]. Based 
on these efforts, it is now accepted that c^S~ and the 
associated high n Rydberg states undergo ultrafast pre- 
dissociation, but the relative importance of tunneling and 
other coupling mechanisms are not clear. Furthermore, 
while the dissociation lifetimes of the u = 1 vibrational 
level have converged to ~ 70/5, the dissociation lifetimes 
of the u = vibrational level still differ by orders of 
magnitude between different studies [T8ti2Ql E2l [23] . 

In contrast to the high n Rydberg states, the elec- 
tron correlation effects in the low n Rydberg states 
(e.g. n = 4,5,6) lead to a scenario in which autoion- 
ization becomes a competing and even prominent decay 
mechanism. Here, we utilize a time-domain pump-probe 
scheme to obtain direct measurements of autoionization 
lifetimes for low n Rydberg states and competing neutral 
predissociation lifetime for the u = level. 

In our experiment, a 1.5 mJ near-infrared (IR) pulse 
with a temporal width of 45 fs was divided into a pump 
and probe path. The pump pulse was focused into a 
hollow-core waveguide filled with Xenon gas to create an 
XUV attosecond pulse train (APT) consisting primarily 
of the 13th and 15th harmonics. The temporal duration 
of the APT envelope was approximately 3 fs. The XUV 
pulse was then focused into an effusive gas jet of molec- 
ular oxygen to create a temporally localized wavepacket 
corresponding to Rydberg nlag states. A probe IR pulse 
with a peak intensity of 4 TW/cm^ was used to ionize 
the neutral excited wavepacket resulting in the forma- 
tion of the continuum molecular ion, {c^Ti~)02 - This 
state served as our final product channel as it leads to 
formation of easily observable 0+ ions. We used a veloc- 
ity map imaging (VMI) detector to measure the kinetic 
energy and yield of 0+ ions along two dissociation path- 
ways as a function of pump-probe time delay. 

To understand the final product channel, we first per- 
formed a non-time-resolved study of the (c^I]~)02" state 
dissociation pathways by tuning the 15th harmonic of the 
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FIG. 2. (a) The relevant photodissociation channels of 
(c^E-)0+ which serve as the final product channel in our 
pump-probe measurements discussed later, (b) The kinetic 
energy spectrum of for direct excitation of (c^S~)0^ with 
the 15t/i harmonic. The color-coded Gaussian fits that corre- 
spond to the dissociation channels in (a) are overlaid on the 
kinetic energy plot. The inset shows the reconstructed VMI 
image of ions. 



XUV APT to be directly on resonance with this ionic 
state. The two vibrational levels (y = 0, 1) of (c^I]~)02" 
predissociate along well-known pathways associated with 
characteristic kinetic energy values [24 . The potential 
energy curves for relevant fragmentation pathways are 
shown in Fig. |2](a). We find three prominent peaks in the 
measured kinetic energy spectrum for 0+ ions, shown in 
Fig. [2] (b) along with the full momentum image for these 
predissociation paths (inset). The most prominent peak 
arises from the (5^11" )02" state which is populated by 
the 13th harmonic of our APT. This dissociation channel 
is not relevant to the present study. 

The outer two peaks observed in the 0+ kinetic en- 
ergy spectrum in Fig. [2](b) originate from the [c^Y,~)02 
state. The peak at 1.9 eV arises from one of two decay 
mechanisms, either tunneling through the potential bar- 
rier or coupling to a repulsive curve, to the final limit 
0+('^5')+0(^D) (L2) [24J. Both vibrational levels have 
been found to participate in this dissociation pathway. 
The outermost peak at 2.9 eV is important in our study 
since it arises purely from the = level which predis- 
sociates to the 0+('^5')+0(^P) (LI) limit [25, 26 . 

Now we turn our attention to the measurement of the 
relaxation lifetime for superexcited Rydberg states. We 
tune our XUV pulse such that the peak of the 15th har- 
monic is resonant with the {ns^nd)ag Rydberg states 
(Fig. [sja)). The nla g{c^T^~) Rydberg state assignments 
for the u = vibrational level obtained from prior syn- 
chrotron work [18 are shown in Fig. |3|b) along with our 
excitation spectrum. The excitation of the u = 1 vibra- 
tional manifold is expected to be factor of three times 
weaker than u = [T9] . 

The temporally localized population in the Rydberg 
excited state decays very rapidly due to the strong elec- 
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ion yield in the c^S^ ^Ll channel is given by [27] 
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FIG. 3. (a) A schematic representation of the relaxation life- 
time measurement. The 5scr^ Rydberg state is populated 
by the XUV APT and decays rapidly due to various mech- 
anisms. We interrupt this decay process by removing the 
Rydberg electron with a time-delayed IR pulse to populate 
the (c^E~)Oj state and produce measurable ions with with 
signature kinetic energy values of 1.9 eV and 2.9 eV. (b) The 
XUV excitation spectrum for our 15th harmonic used in the 
lifetime measurements overlaid with the ns and ndag Ryd- 
berg assignments, (c) The ion kinetic energy spectrum with 
the XUV pulse alone (lower curve, blue) as well as in the si- 
multaneous presence of the XUV + IR probing pulse (upper 
curve, red). 



tron correlation effects and the quasi-bound nature of the 
potential energy curve. We interrupt this decay with the 
IR probe pulse to ionize the molecule, leading to the for- 
mation of {c^T^~ which fragments along the known 
pathways discussed earlier (Fig. [5] (a)). The kinetic en- 
ergy spectrum of 0+ is shown in Fig. [s] (c). The lower 
curve represents the background in photo-ion counts due 
to XUV alone. The upper curve shows the photo-ion 
counts in the simultaneous presence of both the XUV and 
IR fields. Clearly, the ion counts in both the LI and L2 
limits are enhanced greatly by the action of the IR probe 
pulse. The relaxation dynamics of the n/<j^(c^S~) Ryd- 
berg states can thus be obtained from the measurement 
of 0+ ion-yield in both the c^T.- ^L2 and c^T.- ^Ll 
dissociation limit as a function of pump-probe time delay. 
As mentioned previously, the yield in the final product 
channel associated with the LI limit comes entirely from 
the u = vibrational level of the c^S~ state. Since the 
excitation from Rydberg states to the continuum ion-core 
is dominated by Au = transitions, the ion yield mea- 
sured in the LI limit arises exclusively from the u = 
population in the Rydberg states. 

The delay dependence of the ion-yield in the final 
product channel is governed by the competition between 
autoionization and neutral dissociation of superexcited 
states. If we consider the rate equations which govern 
the relaxation of a = Rydberg state, we find that the 



yLl{Po,Ta,To]t) = aPo 



-til/Ta + l/To) ^ 



(1) 

where a is the branching ratio for the c^S^ ^Ll chan- 
nel, Po is the initial population in the u = Rydberg 
vibrational level, is the autoionization lifetime, and Tq 
is the predissociation lifetime of the u = Rydberg level. 

This equation can be understood in the context of the 
multichannel quantum defect theory for Rydberg states 
in molecules [28] where the effective quantum number 
(n*) incorporates the effect of core interactions. Us- 
ing this formalism, it has been established that both 
the oscillator strength for excitation to a nlag Rydberg 
state, fni^ and the autoionization rate from that Ryd- 
berg state scale as (n*)~^. However, in the core- ion 
approximation [18], the predissociation rate is indepen- 
dent of the effective quantum number of the Rydberg 
state. The contrasting n* dependence of autoionization 
and neutral dissociation lifetimes is shown in Fig. [4] (a). 
Using Eq. [l] the variation of ion-yield with time de- 
lay can be discussed in three regimes. For low n* states 
(Region I, Fig. [4] (a)), autoionization of the u = pop- 
ulation from the potential well is the dominant mecha- 
nism, manifesting itself as an exponential decay in the 
c^S~ ^Ll channel ion yield (green curve). For high n* 
states (Region III, Fig. [4] (a)), neutral dissociation is 
the dominant mechanism. However, as the neutral frag- 
ments can be photoionized in the dissociation limit there 
is no net loss of experimental ion counts. The delay- 
dependent ion signal therefore exhibits a flat dc response. 
When both mechanisms have comparable rates (Region 
II, Fig. [4] (a)), the exponential decay terminates into 
an appreciable dc baseline. In general, the amplitude of 
the exponential term relative to the dc baseline in Eq. 
[1] elucidates the competition between autoionization and 
neutral dissociation. This allows us to decouple the two 
decay mechanisms and measure the contribution from 
each process. 

Before we quantitatively deduce the autoionization and 
predissociation lifetimes, we note that Eq. [l] refers to 
ion-yield obtained from a single Rydberg level. From 
Fig. [sj^b), we observe that the Rydberg states Ssa^, 
4:d(jg^ 6s(jg 3,11(1 ^dcF g dJCQ prcdomiuaut ly populated by 
the Ibth harmonic. Based on the similarity between 
the energies and effective quantum numbers, it is con- 
venient to classify these Rydberg states [29] into two 
groups i.e. (55,4(i)<j^ and (65,5(i)<j^. Since the har- 
monic field strength is nearly identical for both groups, 
the oscillator strength for each group determines its ex- 
citation probability. Defining r = (^55 4c^)^/(^6s 5d)^ 
the scaling factor for the oscillator strengths and autoion- 
ization lifetimes between consecutive Rydberg states, we 
replace Pq with rPo and with r^/r in Eq. [l] to 
obtain the yield of the second Rydberg group. The 
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FIG. 4. (a) A sketch showing Ta and To scahng with effec- 
tive quantum number n*. The subplots illustrate the ex- 
pected delay dependence of the ion yield in three regimes 
(green curves), (b) Experimental ion yield measured in the 
c^S~ ^Ll channel. The fit obtained for time delays >70 fs 
using the model discussed in the text is shown as a solid green 
line, (c) Experimental ion yield measured in the c^T^u ^Ll 
channel when the harmonic spectrum is tuned to be reso- 
nant with the high n Rydberg states. In this region, neutral 
dissociation dominates resulting in the expected flat delay- 
dependence of the ion yield, (d) Experimental ion yield mea- 
sured in the c^S^ ^L2 channel. Red curves in all figure pan- 
els depict the Gaussian pump-probe cross-correlation profile 
representing additional direct XUV+IR excitation processes. 



time dependence of total yield can thus be written as 
Ylf^'it) = Ylf\P,,T,,To-,t)+Yll''\rPo,Ta/r,To-,t). 
We invoke Rydberg assignments reported by Hikosaka et 
al. [18] to fix the constant r = 0.515. Lastly, we note 
that while we do not include the effects of the higher 
lying Rydberg states, the population of these states is 
substantially smaller (< 10%) and should not affect our 
analysis. 

The delay dependent ion yield obtained in the 
c^S" ^Ll channel is plotted in Fig. [ij^b). The dashed 
red curve depicts the extent of Gaussian cross-correlation 
between the XUV and IR pulses. In our analysis, we 
concentrate on time delays greater than 70 fs where tem- 
poral overlap effects can be ignored. The fitting function 
Yl^^^\t) has three free parameters (Po^Ta^To)^ where Pq 
is an overall multiplier. While the fitting expression con- 
sists of two exponential terms and a dc term, the am- 
plitudes, decay rates and the dc baseline are all inter- 
related through the and Tq parameters. Thus, we 
have a tightly constrained system and the Yl^^\t) fit 



converges to unique values for and r^. The fit shown 
by the solid green curve in Fig. |4jb) yeilds an autoion- 
ization lifetime of = t^sM = 92 ± 6 fs and a neutral 
dissociation lifetime of Tq = 1100 ± 100 fs. Using the 
scaling factor r = 0.515, the autoionization lifetime of 
the second Rydberg group is given by tqsM = 180 ± 10 
fs. These results show that autoionization is the domi- 
nant relaxation process for low to mid n Rydberg states. 
The same understanding can be obtained through a sim- 
ple comparison of the shape of our experimental curve 
(Fig. |4](b)) to the simulated plots shown in Fig. [4]^ a). 

We can look at the other extreme by tuning our har- 
monic spectrum to be on resonance with the high n Ry- 
dberg states near the ionic continuum. Since grows as 
(n*)^, we expect predissociation to dominate, resulting 
in an ion yield that exhibits no dependence on IR pump- 
probe delay (Region III, Fig. [4^ (b)). This is demon- 
strated experimentally in Fig. [4|(c), thus confirming the 
validity of our analysis. 

A similar measurement of the relaxation lifetimes was 
conducted for the c^S~ ^L2 pathway. The ion-yield 
in this channel originates from both vibrational levels, 
complicating the analysis. However, due to a lower ex- 
citation probability and an ultrafast neutral dissociation 
lifetime (n = 69 fs [SI ^ |23]), the contribution of the 
V = 1 level to the ion yield is significantly smaller than 
the contribution from the u = level. Adopting a simi- 
lar fitting procedure as above, the ion yield shown in Fig. 
|4|d) yields autoionization lifetimes of r^s,4d = 105 ± 4 fs 
and rQs,5d = 204 ± 8 fs. The deviation of these results 
from the c^S~ ^Ll results discussed earlier can be in- 
terpreted in terms of the unaccounted contribution from 
the u = 1 level. 

To the best of our knowledge, there are no prior ex- 
perimental measurements of autoionization lifetimes for 
these superexcited Rydberg states. In a recent theoreti- 
cal effort [221129], lifetimes of r^s = 90.0 fs {tm = 92.7 fs) 
and tqs = 178 fs {r^d = 183 fs) have been reported, which 
are in a good agreement with our experimental observa- 
tions. We note this agreement also justifies the Rydberg 
state groupings that we had introduced in our model. 
The measured value for predissociation lifetime To is also 
in agreement with the experimental lower bounds mea- 
sured by Hikosaka et al. [18 and Padmanabhan et al. 
[20] who argued Tq > 600 fs and Tq > 1 ps respectively. 

In summary, we were able to disentangle and mea- 
sure the competing autoionizaton and neutral dissocia- 
tion lifetimes of the (5s, 4:d)ag and (6s, 5d)ag Rydberg 
states of O2*, formed through the XUV excitation of a 
2sau electron. This represents the first direct measure- 
ment of the ultrafast autoionization and predissociation 
dynamics in a neutral superexcited molecule. This topic 
has been of significant interest to the synchrotron com- 
munity over last three decades. The knowledge of neutral 
dissociation and autoionization lifetimes can be used fur- 
ther to deduce the rates of various spin-orbit interaction 
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terms and tunneling processes. Thus, the table-top HHG 
based XUV time-domain technique used here represents 
a significant step forward in the understanding of non- 
adiabatic relaxation mechanisms. In general, the cou- 
pling between electronic and nuclear degrees of freedom 
forms an interesting problem where direct time domain 
techniques could be very illuminating [30, 31 . An im- 
portant facet of the time-domain experimentation is the 
possibility of implementing real-time control, wherein a 
light pulse can be used to modify or change the course of 
the photo-chemical reaction in the transient phase. 
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